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Abstract 

We report measurements of the magnetic penetration depth A in single crystals of PrRu4Sbi2 
down to 0.1 K. Both A and superfluid density ps exhibit an exponential behavior for T < 0.5Tc, 
with parameters A{0)/kBTc = 1.9 and A(0) = 2900 A. The value of A(0) is consistent with the 
specific-heat jump value of AC/^T^. = 1.87 measured elsewhere, while the value of A(0) is consistent 
with the measured value of the electronic heat-capacity coefficient 7. Our data are consistent with 
PrRu4Sbi2 being a moderate-coupling, fully-gapped superconductor. We suggest experiments to 
study how the nature of the superconducting state evolves with increasing Ru substitution for Os. 
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The recent discovery^'^ of the Heavy Fermion (HF) skutterudite superconductor (SC) 
PrOs4Sbi2 has attracted much interest due to its differences with the other HFSC Measure- 
ments of dc magnetic susceptibihty, specific heat, electrical resistivity and inelastic neutron 
scattering showed that the ninefold degenerate J = 4 Hund's rule multiplet of Pr is split by 
the cubic crystal electric field, such that its ground state is a nonmagnetic doublet, sepa- 
rated from the first excited state Fs by ~ 10 K. Hence its HF behavior, and consequently the 
origin of its superconductivity, might be attributed to the interaction between the electric 
quadrupolar moments of Pr'^+ and the conduction electrons. It is thus a candidate for the 
first superconductor mediated by quadrupolar fluctuations, i.e. by neither electron-phonon 
nor, as with other HFSC, magnetic interactions. 

Surprisingly, replacement of Os by Ru, i.e. in PrRu4Sbi2, yields a superconductor with 
Tc ~ 1.25 K with significantly different properties. From the slope of the upper critical field^ 
near Tc, {-dHc2/dT)T^ = 2.4 kOe/K and using the procedure in Ref. Q, we get the effec- 
tive mass of the heavy electrons m* ^ 20me- This contrasts with the value m*/me=45 for 
PrOs4Sbi2, showing that while PrOs4Sbi2 is clearly a heavy- fermion material, PrRu4Sbi2 is 
at most, a marginal HF. Various experimental results suggest that these two materials have 
different order-parameter symmetry. Firstly, there is an absence of a Hebel-Slichter peak in 
the nuclear quadrupole resonance (NQR) data^ for PrOs4Sbi2, while a distinct coherence 
peak was seer^ in the Sb-NQR 1/Ti data for PrRu4Sbi2. Secondly, the low-temperature 
power-law behavior seen in specific heaiji and penetration depth-, and the angular variation 
of thermal conductivity^!, suggest the presence of nodes in the order parameter of PrOs4Sbi2. 
Specifically, Ref s. and Q reveal the presence of point nodes on the Fermi surface (FS). For 
PrRu4Sbi2, however, an exponential decrease in 1/Ti is observed below the Hebel-Slichter 
peak^ and was fit with an isotropic gap of magnitude A{0)=1.5kBTc, where A(0) is the mag- 
nitude of the zero-temperature superconducting gap. Thirdly, muon spin rotation (/xSR) ex- 
periments on PrOs4Sbi2 reveal the spontaneous appearance of static internal magnetic fields 
below Tc, providing evidence that the superconducting state is a time-reversal-symmetry- 
breaking (TRSB) state^. Such experiments have not been performed on PrRu4Sbi2. 

In this paper, we present high-precision measurements of penetration depth A(T) of 
PrRu4Sbi2 at temperatures down to 0.1 K using the same experimental conditions^ as for 
PrOs4Sbi2. Both A(r) and superfluid density Ps{T) exhibit exponential behavior at low 
temperatures, suggesting the presence of an isotropic superconducting gap on the FS. Data 
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are best fit by the parameters A{Q)/kBTc = 1.9 and A(0) = 2900 A, and thus suggest that 
PrRu4Sbi2 is a moderate-couphng, fully-gapped superconductor. The values of A(0)/A;b7'c 
and A(0) are consistent with values^ derived from the specific-heat jump AC/7Tc=1.87, and 
the linear specific-heat coefficient 7. 

n 

Details of sample growth and characterization are described in Ref. The observation 
of the de Haas-van Alphen (dHvA) effect from the same batch of samples^, and the large 
residual resistivity ratio (RRR ~ 76), reflect the high quality of the samples. Measurements 
were performed utilizing a 21-MHz tunnel diode oscillator— with a noise level of 2 parts in 
10^ and low drift. The magnitude of the ac field was estimated to be less than 40 mOe. The 
cryostat was surrounded by a bilayer Mumetal shield that reduced the dc field to less than 
1 mOe. The sample was mounted, using a small amount of GE varnish, on a single crystal 
sapphire rod. The other end of the rod was thermally connected to the mixing chamber of 
an Oxford Kelvinox 25 dilution refrigerator. The sample temperature is monitored using a 
calibrated RUO2 resistor at low temperatures {T^ase to 1.8 K). 

The deviation AA(T) = A(r) - A(0.1 K) is proportional to the change in resonant 
frequency A/(T) of the oscillator, with the proportionality factor G dependent on sample 
and coil geometries. We determine G for a pure Al single crystal by fitting the Al data 
to extreme nonlocal expressions and then adjust for relative sample dimensions'^. Testing 
this approach on a single crystal of Pb, we found good agreement with conventional BCS 
expressions. The value of G obtained this way has an uncertainty of ±10% because our 
sample, with approximate dimensions 0.7 x 0.5 x 0.25 mm^, has a rectangular, rather than 
square, basal area^. 

Figure n (O) shows AA(r) for PrRu4Sbi2 as a function of temperature in the low- 
temperature region. The inset shows AA(r) for the entire temperature range. The value of 
Tc, taken to be the mid-point of the transition, is 1.25 K. The data points fiatten out below 
0.22 K (~ O.lSTc), implying activated behavior in this temperature range. Contrasting this, 
we superpose data for PrOs4Sbi2 on the same figure (x), which show A A varying strongly 
(~ T^) with temperature, indicative of low-lying excitationstSi. We fit the PrRu4Sbi2 data to 
the BCS low-temperature expression in the clean and local limit, from Tbase (0.1 K ^ 0.08 Tc) 
to 0.65 K 0.5Tc), using the expression 
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with the proportionahty constant and A(0) as parameters. The best fit (sohd hne in Fig. 
is obtained when A(0) = 2.4 K = l.dksTc. This value is larger than the BCS weak-coupling 
value of 1.76/cbTc, suggesting that PrRu4Sbi2 is in the moderate-coupling regime. To check 
the validity of this value of A(0) we make use of the strong-coupling equation o^^i^^ 



T' \ / 



r/AM = 1 + 5.3 (-^1 In(^), (2) 



Ti \ 2 / 



r^c.M = 1 + 1.8(^1 (^ln(^) + 0.5 ) , (3) 



;i + (^m61n(t)-0.26) 
''^^"'"^ l + (£|)2(i.iin(^) + 0.14) ' ^ ^ 

where each t] represents the correction factor to the corresponding weak-coupling BCS value. 
If we take A(0)=1.9fcBTc, then Eq. Q gives the characteristic (equivalent Einstein) fre- 
quency cjo ~ 17 K and Eq. Q gives AC /'yTc=l.9. This value of AC/'yTc agrees excellently 
with the measured value (1.87) in Ref. jsi giving further evidence that PrRu4Sbi2 is indeed 
a moderate-coupling superconductor. 

To extract the superfluid density ps from our data, we need to know A(0). For a type-II 
superconductor, A(0) can be obtained fron>i^ 



where $o=2. 06x10^ G-A^ is the flux quantum, Hc2{0) is the upper critical field 
at T=0, 6sc=A{0)/kBTc, and 7 is the electronic specific heat coefficient. Using 
{-dH^2/dT)T=2A kOe/K from Ref. S and the expression i7e2(0)=0.693(-d//c2/dr)T,Te, 
we obtain ifc2(0)=2.16 kOe. The superconducting coherence length ,^0 can be estimated 
from the relation iJc2(0)=$o/2vr^o; yielding ^0 ~ 400 A. The value 7=59 mJ/mol-K^ is also 
obtained from Ref. 0. We calculate A(0) using two methods: (1) Taking 5sc=l-9, or (2) 
taking 5sc=l-76 along with strong-coupling corrections (Eqs. |21and|3I). Both methods yield 
A(0) ~ 3200 A. This puts PrRu4Sbi2 in the local limit. Furthermore, from dHvA data^, we 
estimate the mean free path ZRil300 A, implying that the sample is close to the clean limit. 
To calculate ps for an isotropic s-wave superconductor in the clean and local limits we use 
the expression 



where / = [exp(E/A;BT)+l] ^ is the Fermi function, and E = [e'^ + A(T)^]^/^ is the quasi- 
particle energy. The temperature-dependence of A(T) can be obtained by usingi^ 



where Ssc is the only variable parameter, Tc = 1.25 K, a = 2/3, and the specific heat jump 
AC/C = AC/'yTc = 1.87 is an experimentally obtained valued. 

Fig. 121 shows the experimental (Q) and calculated (solid and dotted lines) values of ps 
as a function of temperature. The best fit from 0.1 K to 0.95 K (~0.8Tc) is obtained when 
A(0)=2900 A and A{0) = 1.9kBTc (sohd line). This value of A(0) is 10% lower than the 
earlier-calculated value of 3200 A, but is acceptable because of the uncertainty in obtaining 
the calibration factor G. The value of A(0) once again agrees with the specific-heat jump 
AC/'jTc obtained in Ref.jjvia Eqs. (0) and 0, though it disagrees with the weak-coupling 
value of LbhsTc in Ref. jsl The dotted line in Fig. |2l calculated using the weak-coupling 
parameters A(0)=1.76/cbTc and AC/7Tc=1.43, clearly does not fit the data. Our super- 
fluid data once again suggest that PrRu4Sbi2 is a moderate-coupling superconductor with 
a superconducting gap on the entire FS. 

It is apparent in Fig. |21that the data deviate from local BCS expression above 0.95 K. We 
also noticed a small hump in A A near 1.1 K shown in Fig. ^ which shows up as a curvature 
change in ps in Fig.|21 These features may be due to inhomogeneities and/or impurity effects 
in the sample, because the Tc's of these single crystals varj^ between 1 K and 1.3 K. 

It is puzzling that the substitution of Ru for Os (same column in the periodic table) 
causes PrRu4Sbi2 to differ in so many respects from PrOs4Sbi2. Further, evidence for 
gap anisotropy in the latter compound is contradicted by /iSR measurements^^ (suggesting 
either s or p-wave pairing, with A(0)=2.1A;bTc), scanning tunneling spectroscopy^^, and 
NQR measurements^ (with A(0)=2. GfcsTc). If both PrRu4Sbi2 and PrOs4Sbi2 have isotropic 
gaps, then they are unique, especially the latter, among HFSC, suggesting the possibility of 
(a) an important difference in superconducting properties between HFSC with magnetic and 
non-magnetic /-ion ground states, and (b) a correlation between pairing symmetry (isotropic 
or nodal gap) and mechanism (quadrupolar or magnetic fluctuations) of superconductivityi^. 

Recently, Frederick et al. performed X-ray powder diffraction, magnetic susceptibility 
and electrical resistivity measurementsi^ on single crystals of Pr(Osi_i.RuK)4Sbi2. They 
found that (1) the lattice constant a decreases approximately linearly with increasing Ru 
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concentration, (2) the splitting between the ground and first excited state increases mono- 
tonically with x, with the fits consistent with a doublet ground state for all values of 
X, although reasonable fits can be obtained for a Fi ground state for x near and 1, and 
(3) Tc decreases nearly linearly with substituent concentration away from x=0 and x=l, 
but exhibits a deep minimum (0.75 K) at x = 0.6. The smooth evolution of a and 
with X, and the presence of superconductivity for all values of x, may suggest that both 
PrOs4Sbi2 and PrRu4Sbi2 possess the same order-parameter symmetry. The minimum in 
Tc at X = 0.6 could simply mark the shift from quadrupolar- mediated heavy fermion su- 
perconductivity to phonon-mediated BCS superconductivity. On the other hand, one still 
has to contend with measurement3^i^*2S that indicate point-node gap structure, TRSB and 
double superconducting transitions in PrOs4Sbi2. If so, the minimum in could be a con- 
sequence of competing order-parameter symmetries with a possible quantum critical point 
between them. It is also interesting to notice in Ref. |l9| that the step-like structure seen in 
ac susceptibility data in the x=0 sample, i.e. PrOs4Sbi2, indicative of an intrinsic second 
superconducting transition, is not seen for all other values of x. 

To further elucidate the relationship between PrOs4Sbi2 and PrRu4Sbi2, work is under- 
way looking at the changes in penetration depth on Pr(Osi_^.Rua.)4Sbi2 for a range of values 
of X. We want to know at which value of x (paying close attention to the value x = 0.6), if 
any, does the isotropic superconducting gap evolve into a nodal one. Moreover, we want to 
know whether the second superconducting transitioit^ seen in PrOs4Sbi2 can also be seen 
in samples where < x < 1. If we do not see the second superconducting transition in 
the other samples, then it is possible that the superconductivity in the two superconducting 
phases of PrOs4Sbi2 respond differently to impurities (Ru substituents) — one is completely 
destroyed by even tiny amounts of impurities, while the other persists (though weakened) 
all the way to PrRu4Sbi2. It might be useful also if the impurity introduced is of an element 
that would not produce an isostructural superconducting compound. 

In conclusion, we report measurements of the magnetic penetration depth A in single 
crystals of PrRu4Sbi2 down to 0.1 K using a tunnel-diode based, self-inductive technique at 
21 MHz. Both A and ps exhibit an exponential behavior for T < 0.5Tc, with parameters 
A{0) /ksTc = 1.9 and A(0) = 2900 A. The value of A(0) is consistent with the specific-heat 
jump value of AC/'-fTc = 1.87 measured elsewhere^, while the value of A(0) is consistent with 
the measured valued of the electronic heat-capacity coefficient 7, and A(0). Our data are 
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consistent with PrRu4Sbi2 being a moderate-coupling, fully-gapped superconductor. We also 
suggest further experiments that can be done to study how the nature of the superconducting 
state evolves with Ru substitution, to further elucidate the relationship among HF behavior, 
superconducting gap structure, and the presence of TRSB in the superconducting state. 
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FIG. 1: (O) Low-temperature dependence of the penetration depth AA(T) of PrRu4Sbi2. The 
sohd hne is the fit to Eqn. nfrom 0.1 K to 0.65 K, with A(0) = 2.4 K = IMbTc. Inset shows 
AA(r) of PrRu4Sbi2 over the full temperature range, (x) AA(T) data of PrOs4Sbi2, taken from 
Ref. a 
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FIG. 2: (O) Superfluid density ps{T) = [X^{0)/X^{T)] calculated from AA(T) data in Fig.IH Lines: 
Ps{T) calculated from Eqn. IHlwith parameters A(0)/A;bTc = 1.9 and ^Cj^T^ = 1.87 (solid line), 
A{0) /ksTc = 1.76 and AC/jTc = 1.43 (dotted line). 
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